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■-1- ABSTRACT.  The  Rice  distribution  arises  in  bombing  problems  where 
referred  to  as  the  "offset  circle  probability  function."  It 
IflSff occurs  in  radar  as  a  basic  relationship  involving  probability  of 
tjjigje  alarm  (PFA),  probability  of  detection  (P<j),  and  signal-power-to- 
^TBtse-power  (S/N)  ratio.  Evaluation  of  Pj  when  PFA  and  S/N  are  known, 
or  of  S/N  when  PFA  and  P^  are  known,  is  of  usefulness  frequently,  but 
available  means  in  the  form  of  curves,  approximations,  and  tables  are 
not  very  satisfying.  This  report  presents  the  solution  of  the  inte¬ 
gral  in  a  form  that  is  interesting  and  easily  programmable.  Programs 
in  BASIC  and  for  the  Hewlett  Packard  9IOOA  desk  calculator  are  also 
i ncl uded. 
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INTRODUCTION 


The  Rice  distribution  arises  in  bombing  problems  where  it  is  re¬ 
ferred  to  as  the  "offset  circie  probability  function."  It  also  occurs 
in  radar  as  a  basic  relationship  involving  probability  of  false  alarm 
(PFA),  probability  of  detection  (P<j ) ,  and  signa l-power-to-noise-power 
(S/N)  ratio.  Evaluation  of  Pj  when  PFA  and  S/N  are  known,  or  of  S/N 
when  PFA  and  P<j  are  known,  is  of  usefulness  frequently,  but  within  the 
writer’s  knowledge,  available  means,  in  the  form  of  curves,  approxi¬ 
mations,  and  tables  are  not  very  satisfying.  The  solution  of  the 
integral  in  a  form  that  is  interesting  and  easily  programmable  is  pre- 
.  sented  herein.  Programs  in  BASIC  and  for  the  newly  available  Hewlett 
Packard  9I00A  desk  calculator  are  also  included.  t 


DISCUSSION  OF  INTEGRAL  VARIATIONS 

The  integral  that  concerns  us  appears  with  numerous  minor  varia¬ 
tions  in  several  references.  Six  of  these  variations  are  as  follows: 

Equation  I  is  the  offset  circle  probability  function. 


P(a,v) 


[  xe'(a2  +  x2>/2  lQ  (ax)  dx 

"o 


(I) 


Reference  I  on  page  15,  Eq.  1.5  gives 


Reference  2  on  page  407,  substituting  Eq.  7-144  into  Eq.  7-143  gives 


In  Eq.  2 


E  =  r,  dE  =  dr 
n  n 

Et  =  b'  Pd  =  Pd 

ip  =  N,  E  2/2<p  =  S2 
*  s  ro 


In  Eq.  3 


In  Eq.  4 


P  =  P 
sb  rd 


e  =  r,  e*  =  b 


t2  = 


=  n,  j  =  sSm,  d  =  p , 


In  Eq.  5 


V-  '°9e  P 


e  nb 


b 

«£n 


y2  =  r2/2N,  y  =  r/y*§N,  dy  =  dr/,/5N 

In  Eq.  6  (essentially  the  same  as  2) 

R  =  r,  VT  =  b 

ip  =  N,  A2/2<p  =  S2 

o  '  o 


In  Eq.  1-7,  the  modified  Bessel  function  of  the  first  kind  and 
zero  order  l0(Z)  appears,  and  its  expansion  in  series  form,  as  found  in 
handbooks  or  in  Ref.  1-5,  is  given  as  fol.lows: 


00  72n 

I  (Z)  =  E  - - - 

0  n=0  22n(n!)2 


(8) 
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Substituting  Eq.  8  into  Eq.  7  and  replacing  Z  with  rs  /2/N  gives 


P 


d 


(r  s  sm)zn 

22n(n!  )2 


(9) 


In  Eq.  9,  since  we  are  integrating  with  respect  to  r,  and  summing 
with  respect  to  n,  we  can  take  e"s2/N  outside  the  integral  and  take 
re"f2/2N  inside  the  summation  to  give  the  following: 


P 


d 


(s^7N)2n  e"rV2N  r2n+1 


,2n 


(n!)2 


dr 


(10) 


In  Eq.  10,  the  s  and  r  functions  that  appear  inside  the  summation 
sign  have  also  been  separated. 

The  integration  and  summation  in  Eq.  10  can  be  interchanged  in  order 
and,  since  we  are  integrating  with  respect  to  r,  we  can  write  the 
fol lowing. 


In  Eq.  II  consider  the  integral 


2n+l 

r 


dr 


(12) 


Let 


y  =  -r2/2N  or  r2  =  -2Ny 


A 


Then 


or 

dr  *  -  p-  dy 


By 


making  these  substitutions  into  Eq.  12  we  get  Eq.  13. 


a 

*'r=b 


-N  I  ey(-2Ny)n  dy 


Expressing  Eq.  13  operationally  we  get 


^iev(-2Ny,n]  ^ 


r=* 


r=b 


N  ey  (1  -  D  +  D2  -  D3  +  D4  -  D5 - >(-2Ny) 


4 

=  j^-N  ey{(-2Ny)n-n(-: 


*»  ,, 

r=°° 

n 

r=b 


2N)(-2Ny) 


n-l 


+<n)(n-l>(-2N)2(-2Ny> 


-h]:; 


(13) 


(14) 


(15) 


(16) 


(17) 
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The  form  of  £q.  24  is  noteworthy  in  that  for  each  value  of  n,  the 
second  summation  is  the  ratio  of  the  sum  of  the  first  n  terms  of  the 
series  divided  by  the  entire  e^  series,  and  the  term  between  the  summation 
signs  acting  as  a  factor  is  the  ratio  of  the  nth  term  of  the  eR  series 
divided  by  the  entire  eR  series. 

in  programming  Eq.  24,  it  will  help  to  write  out  the  development 
of  for  the  f i rst  few  va 1 ues  of  n . 

n  k  Partial  Sums 

0  0  (l/eR)(l/eT) 

1  0,  I  (R/eR)  {  (I  +  T)/eT} 

2  0,1,2  {R2/(eR2!)}  {(I  +  T  +  T2/2!)/eT} 


n  0,  I, - n  <Rn/(eRn«|  {( I  *+  T  +  Tz/2!  + - Tn/n!)/eT( 

We  are  to  find  P<j  with  the  following  given: 

(S/N)  in  db  =  10  log|Q  R 

PFA  =  i/eT 

We  define  and  use  the  following  calculations  where  n  takes  integral 
values  starting  with  zero. 

FI  =  Rn/(e?n!) 

F0  =  Tn/(eTn!) 

F2  =  (I  +  T  +  T2/2!  +  -  Tn/n! )/eT 

6  =  iFl )(F2) 

Each  of  the  partial  sums  above  is  made  up  of  two  factors:  the 
factor  involving  R  and  the  factor  involving  T.  We  call  the  R  factor 
"Fln  and  the  T  factor  "F2.n  Succeeding  FI  factors  are  obtained  by 
multiplying  the  preceding  FI  factor  by  R/n.  Succeeding  F2  factors  are 
obtained  by  adding  a  new  term  to  the  preceding  F2  factor;  we  call  this 
additive  increment  "F0.n  Succeeding  F0  increments  are  obtained  by 
multiplying  the  preceding  F0  increment  by  T/n. 
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In  the  program,  we  call  the  partial  sums  *'6”  because  they  are 
additive  increments  to  the  previous  sum  to  get  the  next  approximation 
to  P(j.  These  6  terms  ''an  be  very  small  for  initial  values  of  n  and 
therefore  we  cannot  use  a  simple  comparison  of  a  newly  calculated  6 
with  a  small  number  as  a  test  for  stepping  the  calculation.  We  com¬ 
pare  each  6  with  the  previous  6  and  continue  as  long  as  6s  are  in¬ 
creasing.  When  the  6s  start  decreasing  we  terminate  when  we  reach 
a  6  £  I0*1 1  value,  which  goes  beyond  needed  accuracy  and  can  be  altered 
if  desired,  although  the  required  number  of  executions  will  not  be 
greatly  reduced. 

A  program  for  Eq.  24  in  BASIC  is  as  follows: 


10 

READ 

PO,  P.4  Cwhere  P0=PFA,  R4=(S/N)  in  db;  both  given] 

20 

LET 

T=  -L0G(P0)  Cwhere  LOG  is  to  base  e] 

30 

LET 

R=l0t( . I*R4) 

40 

LET 

P4=0  Cwhere  P4=Pcj ;  to  be  calculated] 

50 

LET 

D 1  =0  Cwhere  D i =6 ;  see  02  below] 

60 

LET 

M=0  Cwhere  M=n] 

80 

LET 

FI=EXP(-R) 

100 

LET 

F2=0 

120 

LET 

F0=P0 

160 

LET 

F2=F2  +  F0 

180 

LET 

D2=FI#F2  Cwhere  D2=6,  Dl  and  D2  are  succeeding  values  of  6] 

200 

LET 

P4=P4  +  D2 

220 

LET 

M=M  +  I 

240 

LET 

F 1 =(R/M)#F 1 

260 

LET 

F0=(T/M)#F0 

280 

IF 

D2>DI  THEN  340 

300 

IF 

D2> 1 E-l 1  THEN  340 

320 

GO  TO  380 

340 

LET 

D 1  =D2 

550 

GO  TO  160 

380 

PRINT  "PFA=  P0,  "S/N=  R4,  ”PD=  P4 

400 

GO  to  10 
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500  DATA  -  CEn+er  PFA,  In  pairs] 


999  END 

Finding  (S/N)  =  R  when  PFA  =  l/e^  and  P<j  are  given  can  be  done 
as  fol lows: 


»  Rn  n  Tn-k 

f  (R)  =  2  Z  - - 

n=0  e  n<  R=0  e'[(n-k)!] 


=  P.(R)  -  P. 
d  d 


We  wish  to  find  the  value  of  R  f hat  makes  f(R) 
Taking  Pj  and  T  as  constants  we  get: 


=0  in  Eq .  25 . 


, ,  "  Dn  n  xn-k 

=  E  E  “T^ - Pd(R) 

dR  n=0  eR  Rn!  k=0  eT[(n-k)!]  d 


=  l/eR  l/eTC0+l (T+l )+R(Tz/2!+T/l !+l ) 


+  R2/2!(T3/3!+T2/2!+T/l!+l>  - ] 


-  l/eR  l/eTCl+R(T+l)+R2/2!(T2/2!+T/l!+l) 


+  R3/3! (T3/3!+T2/2!+T/l !  +1)  - ] 


df  (R) 
dR 


i  /  R  i  /  T  fr  .R  T2  ,  R2  T3  .  R3  T4  *1 

-  I/e  l/e  [T  +  ,,  2!  +  2!  3!  3!  4!  J 


If  we  make  an  initial  estimate  of  R  =  Rj,  we  can  use  the  Newton- 
Raphson  iterative  procedure  (Eq.  28)  to  get  succeeding  approximations. 


R2  =  Ri  -  f  (Rj.  )/f 1  (Ri ) 
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The  function  f(R)  plotted  as  the  ordinate  with  R  as  the  abscissa 
lies  between  f(R)  =  PFA  -  P<j  and  an  upper  asymptote  f(R)  =  I  -  P<j; 
at  R  =  0,  f(R)  =  PFA  -  P(j  with  a  slope  T/e^  (where  1/eT  =  PFA)  and  rises, 
changing  from  concave  upward  to  concave  downward,  toward  its  upper 
asymptote.  Because  of  this  shape,  the  repetitive  use  of  Eq.  28  will 
converge  only  when  the  initial  Ri  has  suitable  values.  Two  programs, 
one  in  BASIC  and  one  for  the  H.P.  9I00A,  are  presented  below.  In  the 
BASIC  program,  a  rough  initial  estimate  of  R  is  generated.  In  the  H.P. 
9I00A  program,  a  value  for  R*  must  be  entered  and  changed  if  necessary 
to  secure  convergence.  The  BASIC  program,  as  written,  finds  S/N  in  db 
to  the  nearest  0.0 i  db.  The  H.P.  9I00A  program  presents  the  result 
following  each  iteration  showing  Rj  and  R2  for  comparison;  pressing 
•'CONTINUE”  proceeds  to  the  next  approximation  so  that  the  full  accuracy 
of  the  calculator  is  attainable.  When  the  two  values  are  adequately 
alike,  setting  the  flag  and  pressing  "CONTINUE"  gives  a  display  of  S/N 
in  decibels.  A  separate  program  for  the  H.P.  9I00A  is  given  also  to 
generate  an  initial  value  for  Rj.  The  BASIC  program,  as  given,  has  been 
run  on  the  GE  .Mark  I  time-sharing  system  for  a  wide  range  of  PFA  and  Pq 
values,  but  it  operates  for  only  a  limited  range  of  values  on  the  GE 
Mark  II  time-sharing  system  on  which  it  produces  underflows  that  some¬ 
times  lead  to  divisions  by  zero. 


CALCULATION  OF  S/N  WHEN  PFA  AND  Pd  ARE  GIVEN 

All  alpha  registers  and  register  9  are  used  for  storage.  If  new 
estimate  R2  diverges,  insert  an  intermediate  estimate  in  the  y  register 
before  pressing  "CONTINUE"  at  program  step  71  (not  necessary  to  preserve 
x  and  z  displays).  If  R2  is  negative,  zero  can  be  used  as  a  new  R2. 


ktuin  *rrv  *  aa« 

itno  » r  nyyy 


- 1 - 

UtfUf  1 

Mara*  | 

s*»  j  *t» 

Ml 

X 

y 

X 

f 

* 

d 

e 

b 

•  ! 

«:oj"f 

60 

•i'& 

43 

6;2:« 

06 

6*  3  i7 

07 

6 ;  4  GO  TO 

44 

6!  5  1 

01 

« :  6  > 

13 

6 '  7  RCl 

61 

•[?.]-  1 

34 

<* 

6  *  a  '♦ 

I  I 

33 

... _ 

.  .  j 

•■Aji.  _ 

17 

*\c  ;♦ 

(35 

*:</  4 

25 

1 7 :  o  “  I34 

1  7 :  1  (STOP 

l« 

*z 

*1 

New  f 

stlna 

*  of  (S/k) 

r  i 

.appears 

* 

15 

In  y 

regls 

1  i  r 

1 

«r  with  prl 

! 

ior  esti- 

7;  3  it 

27 

♦ 

1 

■ate 

In  th 

i  z  registei 

Press 

7 '4  :il,T 

;  x. 

;64 

•coin 

JNUE" 

i  | 

to  go  to  next  esti- 

7.5!  - 

!34 

■ate. 

Mien 

i  i  1  i 

Ri  and  R7  are  ade- 

7;6  1 

01 

L  .  . . 

quate 

ly  alike,  aanualiy  set' flag! 

7 ;  7  ♦ 

33 

r 

and  press 

■CONTINUE*  to  convert 

7  8  y*() 

40 

to  de 

i  ;  — 

■cibels:  see  Iasi 

1  '  *  !  • 

t  program 

7-9  f 

15 

;  . 

line. 

i 

7  &  ♦ 

25 

| 

7;6  « 

13 

i 

7.c;£$j„ 

32 

r 

L  . . . 

7  t/'e* 

74 

i 

18 


APPROXIMATION  TO  S/N  WHEN  PROBABILITY  OF  FALSE  ALARM  AND 
PROBABILITY  Oh  DETEC  I  ION  ARE  GIVEN 

This  is  intended  onjy  for  use  as  a  first  estimate  of  R  in  Eq.  28. 
APPROXIMATE  R  =  /T  (A  +  2/2(2Pd  -  I))  +  8  Pd(Pd  -  I)  +  1.5 
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A  program  in  BASIC  for  Eq.  28  is  as  follows: 

10  PRINT  "PFA,"  "PD,"  "(S/N),  DB" 

15  FRINT 

20  READ  PO,  P4 

25  IF  PO  =  P4  THEN  365 

30  IF  PO  >  P4  THEN  378 

35  LET  T  =  -LOG(PO) 

40  LET ,RI  =  T  +  SQR  (8#T)»(2  P4  -i)  +  8*  P4#(P4 

49  LET  R2  =  Rl 

50  IF  Rl  >  0  THEN  70 

60  LET  Rl  =  (P4/P0  -  l)/T 

70  LET  N  =  I 

75  LET  F  =  0 

80  LET  FO  =  PO 

90  LET  FI  =  EXP(-RI) 

100  LET  S  =  0 

I  10  LET  D I  =  0 

120  LET  D  =  F l*P0 

130  LET  P  =  D 

140  LET  FO  =  (T/N)*FO 

150  LET  SO  =  FI*F0 

160  LET  S  =  S  +  SO 

170  LET  D  =  (RI/N)#(S0  +  D) 

180  IF  Dl  <  D  THEN  210 

190  IF  IE-1  I  <  D  THEN  210 
200  LET  F  =  I 

210  LET  Dl  =  D 

220  LET  P  =  P  +  D 

230  LET  FI  =  (Rl /N)*F I 

240  LET  N  =  N  +  I 

250  IF  F  =  I  THEN  265' 

260  GO  TO  140 


22 


265  LET  RO  =  ( P4  -P)/S 
280  IF  ABS(RO)  <  10  THEN  300 
290  LET  RO  =  IO*(RO/ABS(RO) ) 

300  LET  R  =  Rl  +  RO 
310  IF  R  >  0  THEN  330 

320  LET  R!  =  RI/2 

325  GO  TO  70 

330  IF  ABS  (L0G(R/R1 ) )  <  =  .00115  THEN  360 
340  LET  Rl  =  R 

350  GO  TO  70 

360  PRINT  PO,  P4,  INTdOOO  *  (LOG(R)/LOG(  10) )  +  .5)/IOO 

361  PRINT  R2,  R 

362  GO  TO  15 

365  PRINT 

370  PRINT  PO,  P4,  "NO  SIGNAL" 

375  GO  TO  15 

378  PRINT 

380  PRINT  PO,  P4,  "PD  <  PFA" 

385  GO  TO  15 

500  DATA  PFA,  Pd,  -  PFA,  Pd,  ETC. 

999  END 

CONCLUSIONS 

These  programs  have  proven  useful  in  evaluating  radar  data  acquired 
using  pulse  height  analysis. 


NWC  TP  4999 


REFERENCES 


1.  Barton,  David  K.  Radar  System  Analysis.  Clifford  Heights,  N.J., 
Prentice-Hall,  1964. 

2.  Schwartz,  Mischa.  Information  Transmission ,  Modulation, ,  and  Noise. 
New  York,  McGraw-Hi 1 1,  1959. 

3.  Wainstein  L.  A.,  and  V.  D.  Zubakov.  Extraction  of  Signals  from 
Noise.  Clifford  Heights,  N.  J.,  Prentice-Hall,  1962. 

4.  Berkowitz,  Raymond  S.  Modem  Radar:  Analysis,  Evaluation,  and 
System  Design.  New  York,  Wi ley,  1965. 

5.  Skolnik,  Merrill  I.  Introduction  to  Radar  Systems.  New  York, 
McGraw-Hill,  1962. 


24 


pa™™™ 


DOCUMENT  CONTROL  DATA  .  R  &  D 

(Security  c teeelticetien  mi  title,  body  of  mb*  tree  t  end  Indealnd  «*n notation  must  be  entered  when  the  ovetetl  report  Ik  clateilledl 


».  ORIGINATING  ACTIVITY  (Carptete  mother)  Zm,  REPORT  SECURITY  CLASSIFICATION 

Naval  Weapons  Center  UNCLASSIFIED 


China  Lake,  California  93555 


3  REPORT  TITLE 


SOLUTIONS  AND  SIMPLE  PROGRAMS  TO  PERFORM  ACCURATE  CALCULATIONS  BASED  ON 
RICE  DISTRIBUTION 


4.  DESCRIPTIVE  NOTES  (Type  of  report  end  /ncfuiir*  detee) 


S.  au  tmorisi  (Fleet  nee re,  middle  inlttel.  Uetneme) 

Jerome  P.  Madden 


S-  REPORT  DATE 


September  1970 


SO.  CONTRACT  OR  GRANT  NO. 


*.  PROJECT  NO. 


70.  TOTAL  NO.  OP  PACES  7*.  NO.  OP  REPS 

24  5 


•O.  ORIGINATOR'S  REPORT  NUMRCRlSl 

NWC  TP  4999 


AirTask  A30-303/2 16 l/W 1 174-000-01 


10.  DISTRIBUTION  STATEMENT 


EACH  TRANSMITTAL  OF  THIS  DOCUMENT  OUTSIDE  THE  AGENCIES  OF  THE  U.S.  GOVERNMENT  MUST 
HAVE  PRIOR  APPROVAL  OF  THE  NAVAL  WEAPONS  CENTER.  ^  ca  ur  mC  u*s*  GOVERNMENT  MUST 


12.  SPONSORING  MILITARY  ACTIVITY 

Naval  Air  Systems  Command 
Naval  Material  Command 
Washington,  D.C.  20360 


IS.  ABSTRACT 


■  ''The  Rice  distribution  arises  in  bombing  problems  where  it  is 
referred  to  as  the  "offset  circle  probability  function."  It  also 
occurs  in  radar  as  a  basic  relationship  involving  probability  of 
false  alarm  (PFA),  probability  of  detection  (P^),  and  signal- 
power-  to-noise-power  (S/N)  ratio.  Evaluation  of  P£  when  PFA  and 
S/N  are  known,  or  of  S/N  when  PFA  and  P£  are  known,  is  of  usefulness 
frequently,  but  available  means  in  the  form  of  curves,  approxima¬ 
tions,  and  tables  are  not  very  satisfying.  This  report  presents  the 
solution  of  the  integral  in  a  form  that  is  interesting  and  easily 
programmable.  Programs  in  BASIC  and  for  the  Hewlett  Packard  9I00A 
desk  calculator  are  also  included.  ( 


.1473 


(PAGE  1) 


S/N  0101- 807-6801 


urity  cation 


Rice  Distribution 
Modified  Rayleigh  Distribution 
Offset  Circle  Probabiliry  Function 
Solution  of  Integrals 
Probability  of  Detection 
Probability  of  False  Alarm 
Signal-to-Noise  Ratio 
Computer  Programs 


DD  1473  iwcxi 

(PAPE  2) 


